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The Generalized Kohn Theorem in arrays of parabolic quantum dots with Coulomb interaction 
between the dots is violated, if there are different dot species involved. We solve the Schrodinger 
\ equation for cubic lattices with two different dots per unit cell: i) two different circular dots and 

■ ii) two elliptical dots, which are rotated by 90° relative to each other. The interaction between the 

dots is considered in dipole approximation and long- wavelength excitation spectra including FIR 
intensities are calculated. The energy spectrum of the first case can be expressed as a superposition 
of two noninteracting dots with an effective confinement frequency, which includes the effect of dot 
interaction. Only in the second case a splitting of degenerate absorption lines and an anticrossing 
occurs, which is a qualitative indication for interdot interaction. If the interaction becomes very 
strong and if all lattice sites (not necessarily confinement potentials) are equivalent, then the con- 
tribution of the dot interaction outweigh possible differences in the confinement potentials and the 
Kohn Theorem gradually reentries, in the sense that one pair of excitation modes (pseudo Kohn 
modes) becomes independent of the interaction strength. 



a 

C3 , PACS: 73.20.D (Quantum dots), 73.20.Mf (Collective Excitations) 



Ch 1 I. INTRODUCTION 

o 

The Generalized Kohn TheoremEl (GKTh) plays a crucial role in quantum dot physics with far reaching conse- 
quences. It considers interacting electron systems in a harmonic confinement and a constant magnetic field, and it 
states that excitations by long wavelength radiation are not effected by the electron electron (e e) interaction. This 
statement applies to arrays of identical harmonic dot confinements (with e e interaction between the dots) as well (see 
' Appendix of|^). This does not mean that all excitations are independent of e e interaction, but only the optically active 
ones (Kohn modes), and it does not mean that e e interaction is not important for the other excitations. However, this 
C — ■ fact prevents the e e interaction from beeing seen and investigated e.g. by far infrared (FIR) spectroscopy. The FIR 
absorption spectrum of the whole system agrees exactly with the spectrum of a single particle. The GKTh does not 
hold for arrays of different dot confinements, e.g. periodic dot lattices with two different harmonic dot confinements 
per unit cello. Then, all collective modes are excited by FIR radiation and effected by e e interaction, or in other 
words, there is no Kohn mode. The calculation and investigation of absorption frequencies and probabilities in the 
latter case is the subject of this work. 

In order to obtain a visual picture, let us first consider a classical model for the Kohn mode for vanishing magnetic 
• field. (This preliminary consideration will be replaced by a rigorous quantum mechanical treatment in the following.) 
Classically, the charge distributions of all dots oscillate rigidly in- phase with the bare confinement frequency, and 
O ■ the e e interaction contributes only a constant term to the total energy (independent of elongation) . If we have more 
than one identical dots per unit cell, there are additional collective modes, in which the individual dots oscillate out of 
phase, and which are affected by dot interaction, but which are not optically active. Consequently, the dot interaction 
is not observable with FIR spectroscopy in arrangements of identical dots. One way to trick Kohn's theorem is to 
include different dot species. Then, there is no coherent oscillation mode for all dots, which docs not change the 
■ e e interaction energy of the system in elongation, because there is no common bare confinement frequency. As a 
consequence, all collective modes (two modes per dot in the unit cell) are effected by dot interaction and excited by 
FIR radiation with a finite probability. In other words, the Generalized Kohn Theorem for dot arrays is broken. 
Other systems, where Kohn's Theorem does not hold, comprise: i) anharmonic confinementsu'u (circular dots with r 4 
and higher order terms in the radial dependence or cubic dots with terms of type x 2 y 2 ) , ii) hole dots with different 
effective massesll One point of this paper is that the GKTh can be broken despite an exactly harmonic. Hamiltonian. 
A further possibility to observe the e e interaction in the excitations is to consider finite wave lengthau. 
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II. MAGNETOPHONON HAMILTONIAN 



The first part of the calculation of the eigenstates of the Hamiltonian follows closely the procedure described in Ref. 
||. We only have to consider that now the confinement potentials and electron numbers can be different in different dots. 
After introducing center- of - mass (cm.) and relative coordinates in each dot and applying the dipole approximation 
for the Coulomb interaction between the dots, we observe that the Hamiltonian of all cm. coordinates is decoupled 
from individual dot Hamiltonians in the relative coordinates. That's why all excitations can be classified into i) 
collective (cm.) excitations, and ii) intra-dot excitations. The latter are not considered here because they are not 
optically active. The Hamiltonian in the cm. coordinates R„. Q reads in atomic units ft = m = e = 1 (see also Sect. 
IV A in Ref. g) 



2^ V^ + ^~ A(U ^ ) 

n, a L v 



(i) 



where V n ,a = Rn.a - Rn°a is the elongation of the cm. at lattice site (n, a) and P n . Q = —i Vu„ „ is the corresponding 
canonical momentum operator, n runs over the unit cells and a over the dot species within a cell. N a is the number 
of electrons in dot a, and m* the effective mass. It is clear already from inspection of (0) that the eigenvalues of H cm _ 
do not depend on the explicitly shown electron numbers N a , because the factors y/N a can be considered just as a 
rescaling factor of the coordinates U„. Q . However, the eigenfunctions (and quantities derived from them) do depend 
on the explicit N a . The force constant tensor reads 



a > = -e-WN a N a , T (r(°) - R^ Q ,) for (n, a) ^ (n' , a') 



(2) 



(3) 



where e~ 1 is the inverse background dielectric constant and Q a the bare confinement tensor, which produces a 
harmonic confinement. The dipole tensor is defined as T(a) = \ [3 a o a — a 2 I] where (o) denotes the dyad 
product and I the unit tensor. Observe that C depends on N a implicitely which effects the energy eigenvalues. 
A unitary transformation to collective magnetophonon coordinates 
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where N c is the number of unit cells, leaves us with a sum on N c decoupled subsystems H c 

2m* ^ + - A(U <^ 



Eq H <i 
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A(U* ) 



(6) 



which includes the dynamical matrix 



With (g) and (|), we obtain 



(7) 
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C q;Q , Q = ^«+e 1 N a ^2 T(a a -a a >) 



E 



^T(R(«) + a o -a a ,)-e^ 0, T(Rf; 

_ a.' 

C q;ce , Q - = -e- l ^fN a N a , e Jq ' R " 'T (r[, 0) + a a - a a ,) for a ? a' 



(8) 

(9) 
(10) 



where R< 0) Q = R^ + a Q and n^O under the sum means that the term R^ ' = is excluded. 




©--© 



FIG. 1. Minimum unit cells for the two dot architectures considered in this paper with two different circular dots (left) and 
two identical, but rotated, ellipsoidal dots (right). 



Now, we focus our attention to long- wavelength modes (the index q = is dropped henceforth) and consider a 
simple cubic lattice, alternatively occupied by two different dot species. The minimum unit cell is face centered cubic 
(see Fig.l) with lattice constant a. After performing the lattice sum involved in (^) numerically, we obtain the four 
2x2 dynamical matrices 



Cn = O x + d Pl I 
C 2 2 = J~&2 + dp 2 I 
C12 = C 2 i = -dpi2 I 

with the interaction parameters 

Pi = 2N i e^ 1 /{n.n.distance) 3 = Av^N^ 1 /a 3 , (* = h 2) 
Pl2 = 2y/NiN 2 e~ 1 /(n.n.distance) 3 = Ay/2 y/NiN 2 e~ 1 /a 3 

and d = 1.460. From the preceding definitions it follows that pi 2 = JpxPi- 



(11) 
(12) 
(13) 



(14) 
(15) 



III. EIGENSTATES 



Now we are going to find eigenvalues and eigenfunction of (g). For avoiding divergences for B = 0, we add an 
isotropic oscillator potential \ J2 a ^o^a *° l ne kinetic energy in (^J) and subtract it from the interaction term, luq is 
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in principle arbitrary, but we chose the mean value of the bare confinement frequencies included in fij and f^- Now 
we replace the coordinates in (pi) (for q = 0) by Boson ladder operators. This is analogous to the usual text book 
transformation (see e.g. Ref. gSect. 3.3) apart from the factors y/N a . It is obvious that this modification can be 
taken into account by introducing scaled coordinates U Q — ► U Q = y/N a U Q (what implies P a — > P Q = ~P a /\JN a ). 




O-al + Oa2 



N a U ay — — — a — a al — a a2 — a a i + a a2 



(16) 
(17) 



where the first subscript (a = 1, 2) indicates the dot number and the second one the component. The transformation 
of the cm. momentum operators is analogous. 




a a2 ~ a al — a a2 



a a i — a a 2 



(18) 
(19) 



The cyclotron frequency is lo* = B/m*c and u>* 



'ui* 2 + 4u>q. Firstly, it is clear that the Hamiltonian in these 



ladder operators does not show an explicit N a — dependence anymore (apart from that implicit in the dynamical 
matrix). This implies that the eigenvalues of the Hamiltonian do not depend on those N a explicitly seen in (|^). 
Secondly, the commutators of the ladder operators are not influenced by the iV Q -factors and agree with those of 



Bosons: [a a i, 



,+ 1 



I and all other commutators vanish. (This is because the commutators of the U Q and P a i agree 



with the commutators of the untildcd quantities.) Now, the total Hamiltonian can be written in matrix notation in 
the following compact form 



H = [a+a] ■ H 



a 

a+ 



where 



[a + a] = [af t af 2 aj 2 | an a 12 a 2 i a 22 ] 



(20) 



(21) 



and 



a 



is the transposed and Hermitian conjugate of (VW- The 8x8 Hamiltonian matrix is not unique, but can 



be cast into the following form 



consisting of the 4x4 matrices 



H 



a (i 
(3* a* 



with a + 



0T = f3 



(22) 
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(23) 
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and the 2x2 matrices 

Cfefe = Cfefc — - cjg I 



Cn C12 
C21 C22 



u+ 

w_ 



E* 



(24) 
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(25) 
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(26) 



Finding the eigenstates of the Boson Hamiltonian (EG) is provided by mathematical physics and described in Ref. 



a 



which preserves Boson commutators 



in full detail. The goal is to find a linear transformation 
and diagonalizes H. We shall only summarize the recipe here. 
The eigenvalues are given by E ni ^ n2 ^ ri3 , ni = 53fc ( n k + 5 ) w fe with n k being non-negative integers and u>k — 2 7^ 



with 7^ being the four positive eigenvalues of the matrix H ■ J. The 8x8 matrix J 

unit matrices. All eigenvalues of H • J come in pairs (7^, — 7fc). 
The eigenf unctions of H are constructed as usual for Bosons 



I 



is made up of 4 x 4 



(1-4) 

\ni,n2,n a ,n 4 >= 



'n k \ 



|0> 



(27) 



The four eigenvectors belonging to the positive eigenvalues are written in the form xj. 



of A + are given by the vectors x^, and by the vectors x^ 



The column vectors 



, which are the eigenvectors belonging to — 7/.. 



The eigenvectors have to be properly orthonormalized ■ J ■ xj, = Si, k - Without degeneracy, the orthogonality is 
guaranteed automatically. The inverse of this particular transformation is obtained from A -1 = J • A + • J which 
shows that the linear transformation is not unitary (but unitary in a non- Euklidian metric). 



IV. OSCILLATOR STRENGTH 



Optical oscillator strength between the states \n >= \rii, n%, 713, 714 > and \n' >= |n' 1; n' 2 , n' 3 , 7i 4 > for polarization 
in ?y = (x or y) direction are defined as 

fn,n'; V — 2 772* oj n>n > \ < n\U mtot \n' > | 2 (28) 

where uj n ,n' is the corresponding excitation energy, and \J v -tot is the 77-component of the total cm. of the electrons 
in a unit cell (apart from a constant term). In formulae, this means Utot = + j^-U2, where iV tot = iVi + 

After expressing the vectors U by ladder operators b kl and using (p7j), we obtain the usual selection rules, i.e., only 
one quantum with energy u>k can be absorped or emitted, so that we obtain only four absorption lines. The result for 
the oscillator strength for the four possible transitions (k = 1...4) and for 77- polarization reads 

fk = | S f . ( K + 1) for / absorption 

Ntotti* \ nk \ emission v ; 

where n k denotes the initial state, and 




In the last definition, Uki and Vk% for i= 1 . . .4 are the components of the vectors and Vfc, respectively. The oscillator 
strength defined in (|28| ) fulfill the following exact /-sum rule ^2 k fk, v = jq^- It is worth pointing out that for equal 
electron numbers in either dot (Ni = N 2 = N), the oscillator strength depends explicitly on N (contrary to the 
optical excitation energies). In all figures presented below the oscillator strength are for iVj = N 2 = N. 
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V. RESULTS 



Now the two simplest cases are discussed in more detail: two different circles and two identical, but rotated ellipses. 
The ratio of the two bare confinement frequencies involved in either case is 1:1.5 which means, that the two confinemnet 
frequencies in units of the mean frequency wo are 1.2 lnd 0.8. In our figures, all frequencies (energies) are given in 
units of the mean confinement frequency ujq and the interaction parameters p in units of u>q . The magnetic field is 
given in terms of the effective cyclotron frequency u>* in units of ujq (upper scale) and in Tesla (lower scale). The 
conversion between both scales is provided by 



0.9134 -lO- 2 ^ . , 
m* u)o\a.u.*\ 



(32) 



In our figures we used ojq = 0.2 a.u* = 2.53 meV and m* of GaAs for this conversion. (We want to stress that 
this parameter choice effects only the magnetic field scale and not the curves.) The definitions of the interaction 
parameters (14) for GaAs in more convenient units reads 



2.26 • 10 7 N, 



(n.n.distance[A]) 3 (LdQ[meV\) 
(For a more detailed discussion of order- of - magnitude estimates see Ref. |.) 



(33) 



For two different circular dots with bare confinement frequencies u>i and u>2 and iVi = N2, the absorption spectrum 
and the oscillator strength are shown in Fig. 2. Although all absorption lines are effected by the dot interaction 
(represented by the interaction parameter p), and all modes are optically active, there is no qualitative effect of 
interaction in the position of the absorption lines. The reason can be understood easily. In this particular case, the 
four eigenmodes can be calculated analytically providing 



where 



e//,1.2 



^1,2, 3,4 



1*)). , 



1,2) 



d± 



(wf + wf) (P1+P2) 



(wf p 2 d + w\ Pi d + uj\ wf) 



(34) 



(35) 



(The upper and lower sign belongs to w e //,i and w e ff,2, respectively). Consequently, if we had to interpret an ex- 
perimental spectrum, we could do this using the formula ( |34| ) for non- interacting dots, but with the effective (i.e. 
interaction affected) confinement parameters defined in (|35|). Only if we take the intensities into account, we see some 
qualitative effect. Whereas for non- interacting dots (with p = 0) and for B — the oscillator strength of all modes 
agree (for a single oscillator, / is independent of the oscillator frequency), there is a large difference for interacting 
dots at p = 0.5. This large difference can be understood as follows. In the limit p — > 00, the upper pair of modes 
develops into the spurious Brillouin zone boundary mode, which has vanishing oscillator strength and the sum rule 
has to be fulfilled only by the lower pair (see also the discussion below). 



In Fig. 2 both dot species bare the same number of electrons. Therefore, only one interaction parameter p is 
involved. Calculations with different Ni (and pi) do not show any qualitative difference. In the limit of large p (and 
equal electron numbers) we obtain from (pq) 



J eff,l,2 



2pd 




± 



; pd 



0{p~ 



(36) 



Consequently, the square of the smaller effective confinement frequency (which is the only one giving rise to modes 
with a finite oscillator strength for large p) approaches the mean value of both squared bare confinement frequencies, 
whereas the larger one grows continously for large p. 

In Fig. 3 and 4b we show the results for two identical, but mutually rotated, elliptical dots. Without dot interaction 
(p = 0), we have two doubly degenerate lines. With increasing interaction strength, we observe a splitting of degenerate 
modes and an anti-crossing behavior for finite B. As in the case of circular dots, the oscillator strength at B = for 
non- interacting dots (p = 0) agree for all four modes. The dot interaction lifts this degeneracy. Additionally, we 
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observe at p = 0.5 that the oscillator strength in the limits of small and large magnetic fields is considerable only for 
two of the modes, except in the gap region, where three modes contribute. By comparison of Fig.s 3a and 4b we see 
that the magnetic field for minimum gap (between the second and third mode) increases with increasing p, whereas 
the gap width decreases. Consequently, the location and width of the gap provides information on the interaction 
strength. 

By comparison of Fig.s 2 and 3 with Fig. 4, and more clearly by consideration of formula (|3^ ) and Fig. 5, it becomes 
clear that in either case the lower pair of degenerate modes at B = converges to a constant (the mean square 
bare confinement frequency yj (tuf + uj%)/2, which amounts to 1.02 luq in our numerical example). Even for finite B, 
there are two branches, which converge to a finite (independent) value for p — > oo, or in other words, which become 
independent of p in this limit. At first sight this looks surprising because the e e interaction does not show any 
saturation, if we increase the interaction parameter, but it continues to compress the dot state. However, there is a 
simple visual explanation for this feature: Generally, the dot interaction adds an additional second order contribution 
to the confinement, which has the same symmetry as the lattice, i.e. it is circular for a cubic lattice. For large p, this 
additional term outweighs the bare confinement, and the effective confinement in both dots becomes isotropic and 
equal. Thus, we approach the case of a lattice of identical dots, for which a pair of Kohn modes exists. Because these 
Kohn modes do not exactly agree with the modes of noninteracting dots, we call them pseudo Kohn modes. In a 
sense, the Generalized Kohn Theorem reentries for dot lattices with strong interdot interaction. In other lattices with 
lower symmetry, the effective confinement in the strong interaction limit might be elliptical, leading pseudo Kohn 
modes with a gap at B = 0. The other pair of modes (which diverge for p — ► oo) turns into the in- folded modes at 
the Brillouin zone corner (because the units cell halves if all dots become equivalent). These modes become spurious 
in the long wavelength and the large-p limit and the oscillator strength of them converge to zero. 
In Fig.s 3a and 4b we observe an additional qualitative effect of dot interaction. For isolated elliptical dots we expect 
a gap between the two excitation branches at B = 0. However, for larger p only the pseudo Kohn mode might be 
observable, because the oscillator strength of the BZ boundary mode decrease rapidly. On the other hand, the two 
lower modes for finite p develope out of the degenerate lower mode for p = 0, whereby the degeneracy at B = 
survives. Therefore, at B = it looks as if we had a circular dot. The closing of the gap between the two most 
intensive branches at B — is not a gradual effect proceeding with increasing p, but initiated by symmetry. (For a 
deeper understanding see also the additional figures in Ref. ?.) 
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FIG. 2. Excitation modes (a) and oscillator strength (multiplied with N) for p — (b) and p — 0.5 (c) for a lattice with 
two different circular dots as described in the text. The radius of the circles in (a) is proportional to the oscillator strength and 
provides a rough overview. 
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FIG. 3. Excitation modes (a) and oscillator strength (multiplied with N) for p = (b) and p = 0.5 (c) for a lattice with two 
identical, but rotated elliptical dots as described in the text and shown in Fig.l. The radius of the circles in (a) is proportional 
to the oscillator strength and provides a rough overview. 



9 




1 1 1 1 ' 1 ' 1 ' 1 1 

1 2 3 4 5 

B [Tesla] 

FIG. 4. Excitation modes for a lattice with two different circular dots (a) and two rotated elliptical dots (b) for a large 
interaction parameter (p = 2). The radius of the circles is proportional to the corresponding oscillator strength. 
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FIG. 5. Pseudo Kohn mode at B — as a function of interaction parameter p for a lattice with two perpendicular elliptical 



We have shown that breaking the GKTh by constructing quantum dot lattices with at least two different dot 
confinements per unit cell has experimentally observable consequences. Generally speaking, there are no Kohn modes, 
i.e. interaction independent modes, anymore. In both of the considered cases, the degeneracy in the FIR intensities 
at B — between the upper and lower absorption lines is lifted due to dot interaction. For two mutually rotated 
elliptical dots (per cell), we observe also a splitting of formerly degenerate absorption frequencies and the appearance 
of an anticrossing. For two different circular dots no qualitative effect of e e interaction in the absorption frequencies 
is observed. Instead, the absorption spectrum can be mimiced by two noninteracting dots with modified (effective) 
confinements. We also pointed out that an extensively strong interaction destroys the effect of interaction by producing 
pseudo- Kohn modes. Although this limit cannot be reached experimentally, it might be important to take this 
tendency into consideration. 

Only in the case of two circular dots there is a simple analytical closed form solution. However, with the formulae 
presented above, the absorption frequencies and oscillator strength for any cubic lattice with two harmonic dot species 
can be easily calculated. The only numerical task is to find the eigenvalues of an explicitly given non- Hcrmitian 8x8 
matrix and to perform a special sum over the eigenvector components. 



I am indebted to D.Hcitmann, J.Kotthaus, and H.Eschrig, and their groups, as well as G.Paasch for helpful 
discussion. 



dots per unit cell. 
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